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Based on the Born-Oppemheimer approximation, we divide total electron Hamiltonian in a spin- 
orbit coupled system into slow orbital motion and fast interband transition process. We find that the 
fast motion induces a gauge field on slow orbital motion, perpendicular to electron momentum, in¬ 
ducing a topological phase. From this general designing principle, we present a theory for generating 
artificial gauge field and topological phase in a conventional two-dimensional electron gas embedded 
in parabolically graded GaAs/Ina;Gai_ 3 ,As/GaAs quantum wells with antidot lattices. By tuning 
the etching depth and period of antidot lattices, the band folding caused by superimposed poten¬ 
tial leads to formation of minibands and band inversions between the neighboring subbands. The 
intersubband spin-orbit interaction opens considerably large nontrivial minigaps and leads to many 
pairs of helical edge states in these gaps. 


PACS numbers: 71.70.Ej, 75.76.+j, 72.25.Mk 

I. INTRODUCTION 

Exploring of various topological quantum states is al¬ 
ways one of the central issue of condensed matter physics 
[H-Q. Topological insulators (TIs) Q, a new class of 
solids, posses unique properties such as robust gapless 
helical ed ge o r surface states and exotic topological exci¬ 
tations . The helical edge states of two-dimensional 
(2D) TIs are protected strictly against elastic backscat- 
tering from nonmagnetic impurities. This feature leads 
to dissipationless conducting channels and therefore is 
promising for possible applications in spintronics, quan¬ 
tum information, thermoelectric transport and on-chip 
interconnection in integrated circuit. These novel appli¬ 
cations require large nontrivial gaps, which suppress the 
coupling between the edge and bulk states, leading to 
dissipationless edge transport. For this purpose, there 
is an ongoing search for feasible realizations of various 
narrow gap materials containing heavy elements, e.g., 
CdTe/HgTe/CdTe quantum wells (QWs) 0-0, and Tin 
film [22|. However, fabrication of high-quality samples 
of these proposed structures still remains a challenging 
task, requiring precise control for material growth. 

In this work, we demonstrate that conventional semi¬ 
conductor GaAs/In 2 ;Gai_a;As/GaAs two-dimensional 
electron gas (2DEG) with antidot lattices can be driven 
into the TI phase. The 2DEGs provide a promising play¬ 
ground for realizing TI states with quite large nontriv¬ 
ial gap 20 meV) operating at liquid nitrigen tem¬ 
perature regime, instead of searching new materials con¬ 
taining heavy atoms. We first present a general analysis 
for generating an artificial gauge field in a semiconduc¬ 
tor 2DEG, then we demonstrate band inversion between 
neighboring subbands utilizing antidot lattices created 
by well-developed semiconductor etching technique, and 
generate the TI phase with many pairs of helical edge 
states. This suggests a completely new method to gen¬ 
erate topological phase in conventional semiconductor 


2DEGs without strong spin-orbit interaction (SOI), at 
liquid nitrigen temperature regime. 


II. GAUGE FIELD FROM 
BORN-OPPENHEIMER APPROXIMATION 

First we discuss the emergence of an artificial gauge 
field in a system of electrons in a 2D system described by 
a low-energy single-particle Hamiltonian H = e(k)(To ® 
To + HXj=i ® Tj, where ai and n {i = 1,2,3) 

are Pauli matrices describing the electron spin and the 
conduction (t^ = +1) and valence (r^ = — 1) bands, re¬ 
spectively, and (To = 70 = / 2 x 2 are identity matrices. 
Taking dzi = Akx, do 2 = Aky, dos = M — Bk"^ and other 
dij = 0, we obtain the Bernevig-Hughes-Zhang (BHZ) 
Hamiltonian for 2D TIs 0. Neglecting the band index 
T, and taking di = —aky — jdkx, d2 = akx + jdky, d^ = 0, 
we get the Hamiltonian for a 2DEG with Rashba and 
Dresselhaus SOIs, where a and /3 are the strengths of 
Rashba and Dresselhaus SOIs, respectively. Next, we 
divide the total Hamiltonian into the intra-band, slow 
part iTorb = e(k)cro O tq and the inter-band, fast part 
Hib = E*,j=i dij (k)(Ti which usually arises from the 
SOIs in real materials. The eigenstate of the total Hamil¬ 
tonian can be decomposed into the fast and slow com¬ 
ponents: |T'(k)) = E„</'n(k)|xn(k)), where {|xn(k))} 
are eigenstates of the fast part Hjb and {(()n(k)} de¬ 
scribe the slow part. The fast spin dynamics compared 
with the slow orbital motion allow us to make the Born- 
Oppenhenmer approximation, i.e., neglecting the cou¬ 
pling between different |xn(k)), and derive an effective 
Hamiltonian governing the slow orbital motion (/)„ (k) (see 
Appendix A): 

HnO^) = (Xn(k) |iTorb(k,f) -t-7TlB(k,CT)|Xn(k)) 

= Horb(k,f - A„)-t en(k), (1) 
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where e„(k) = (xn(k) l-ffml Xn(k)) acts as an effective 
potential that seperates different bands {|Xn(k))} and 
An = -i(xn(k) |5k|x«(k)) is a gauge potential in the 
momentum space of the slow orbital motion, due to the 
interband coupling to the fast spin dynamics [H, [^ . 
For the BHZ Hamiltonian, the gauge potential A„ leads 
to an effective Lorentz force F'a;y,n(k) in the momentum 
space perpendicular to the electric field E; 

-Fxy,n(k) = , y ] = Ay,n] 

A\M + Be) 

” ”2[H2fc2 + (M-H/c2)2]3/2’ ^ 

where Sn = ±1 denotes spin up or down state while 
tn = ±1 denotes the conduction or valence band, re¬ 
spectively (n = 1,2,3,4). The Chern number Cn = 
(l/27r)/g2 d^A:F’ 2 ;j^,n = s„tn [1 + sign (M/H)]/2 is ob¬ 
tained by integrating the field strength F^y^n in the Bril- 
louin zone. The sign change in M would induce a change 
of the Chern number by 1, which corresponds to the topo¬ 
logical phase transition . 

For a 2DEG with Rashba and Dresselhaus SOIs, we 
find ^^^(k) = 0, which means that the Chern number 
vanishes in 2DEG with SOIs. Comparing the Hamilto¬ 
nian of 2D TIs to that of 2DEGs with SOIs, one can see 
clearly that 2D TIs posses an additional degree of free¬ 
dom: the band index r. In order to generate the gauge 
field and realize TI phases in a 2DEG, one needs to create 
minibands and band inversion in 2DEGs. 


III. EFFECTIVE HAMILTONIAN AND 
QUANTUM SPIN HALL STATES 

Based on the above designing principle, we will create 
topological phase in conventional semiconductor 2DEG. 
This is the first demonstration of the formation of a TI 
phase in the s-like band systems, i.e., a 2DEG with 
nanostructured antidot lattice shown schematically in 
Fig. 1(a). Nanostructured antidot lattices, consisting of 
periodically arranged holes that are etched in a 2DEG, 
form a strongly repulsive egg-carton-like periodic poten¬ 
tial in a 2DEG [32l - [3^ . This artificial crystals lead to 
a wide variety of phenomena, for instance, Weiss os¬ 
cillation, chaotic dynamics of electrons, the formation 
of an electronic miniband structure and massless Dirac 
fermions. At low temperatures, the mean free path of 
electrons is much longer than the period of antidot lat¬ 
tices ranging from 10 to 100 nanometers. The modu¬ 
lated periodic potential can also be created by electron 
beam lithography electrodeposition and periodic arrays 
of metallic nanodots can be realized on semiconductor 
surfaces. Due to elastic strains producing these dots, a 
sufficiently strong piezoelectric potential modulation re¬ 
sults in miniband effects in the underlying 2DEG [^.1^. 
Very recently, a honeycomb lattice of coronene molecules 
was created by using a cryogenic scanning tunneling 
microscope on a Gu(lll) surface to construct artificial 


graphene-like lattice with the lattice constant approach¬ 
ing 5nm [ 3 ^. 



FIG. 1: (color online) Schematic of the proposed struc¬ 
ture and its energy bands, (a) A GaAs/In^Gai-a,As/GaAs 
parabolically graded QW with an antidot lattice, which can 
be created by etching technique, (b) Band profile and the 
first and second subbands of the parabolically graded QW. 
(c) Brillouin zone folding induced by a triangular antidot lat¬ 
tice. The numbers 1 — 5 denote the first to the fifth Brillouin 
zones of the antidot lattice, (d) Minibands of the antidot lat¬ 
tice from folding the first and second subbands of the QW 
{Q = 2n/a and a is the antidot lattice constant). Note that 
band inversion occurs between neighboring minibands. 

We consider the 2DEG in a GaAs/In^^Gai As/GaAs 
PQW, which was fabricated successfully before [40l - l4^ . 
with a triangular antidot lattice (see Fig. 1(a)). Before 
going to the numerical calculation, we first give a clear 
physical picture for the emergence of a TI phase in this 
2DEG system upon nanostructuring with antidot lattice. 
The simplest description of the 2DEG system is obtained 
by reducing the eight-band Kane model to the lowest 
conduction subbands of the QW (see Appendix B). This 
gives the Hamiltonian for the 2DEG with periodic anti¬ 
dot lattice potential V{x,y)-. 


H = h'^k'^/2m + Ai2Tz + r]Tx{kxa'^ + kya'^) + V{x,y), (3) 

where Ti (i = x, y, z) are Pauli matrices describing the 
first and second QW subbands {Xn{z)} (n = 1,2) of ef¬ 
fective mass m, and cr* {i = x, y, z) refer to the electron 
spin. The second term A 12 T 2 comes from the energy dif¬ 
ference 2 Ai 2 between the first and second subbands [see 
Fig. 1(b)]. The third term r]Tj;{kxCr^ + kya^) describes 
the inter-subbands SOI (ISOI) obtained from the eight- 
band Kane model using the Ldwdin perturbation theory 
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(see Appendix B). The coupling strength ry is 


V=^ {X2{z)\ 



dzE^jz) 

Ef{z) 


+ P{z) 


d,Piz) \ 

E^{z) ) 


\Xi{z)) 


( 4 ) 

where Eg{z) and Ao( 2 :) are the band gap and spin split- 
off splitting in the QW region, Egi{z) = Eg{z) + Ao( 2 :), 
and P{z) is the Kane matrix element. The SOIs in 
2DEGs usually come from the asymmetry of the QWs, 
i.e., Rashba SOI. Surprisingly, the ISOI can appear in a 
symmetric PQW, behaving like a hidden SOI. From Eq. 

one can see that the ISOI arises from the spatial vari¬ 
ations of the bandgap Eg{z), the Kane matrix P{z), and 
the intrinsic SOI Ao(z), i.e., the variation of the concen¬ 
tration of In component, which behaves like an effective 
local electric field. This local electric field would not push 
the electron and the hole states to the left and right sides 
of the QW, but it can induce a considerably large ISOI 
hidden in symmetric QWs. The initial \xi{z)) and final 
states \x 2 {z)) are neighboring subbands having opposite 
parity, while the variations of Eg{z), P{z) and Ao(z) in 
a symmetric QW are odd. This means that the ISOI can 
exist in symmetric QWs. 

The first and second subbands both form minibands 
due to the Brillouin zone folding caused by the an¬ 
tidot lattice. The band inversion could occur be¬ 
tween the two adjacent minibands of the first sub¬ 
band ■i/'-i- = Xi{z)uk,ni(x,y) and the second sub¬ 
band V'- = X 2 {z)uk,n-z{^,y) (see Fig. 1(d)), where 
Uk,n(x, y) is the nth miniband formed by the antidot lat¬ 
tice. We model the triangular antidot lattice potential 
V (x, y) by a periodic potential Vb[l -I- X]i=i cos (Qi • r)] 
with potential height Vq, Qi = (27r/a)(l, 0), Q 2 = 
(27r/a)(-l/2,-V3/2), Q 3 = (27r/a)(-l/2, v%/2) and a 
is the triangular antidot lattice constant (see Appendix 

C). ^ 

To describe the four minibands (two spin-degenerate 
minibands) |i/;+,t), \tp-A) involved in 

the band inversion, we treat other electron and hole 
minibands by Lowdin perturbation theory and reduce 
the eight-band Kane k-p model to the following effective 
Hamiltonian within the basis (|'!/'-i-, t); IV'-i)-)) IV'-i-;4')) 

IV'-,t)): 


M - Bp Ak_ 0 0 

Ak+ -M + BP 0 0 

0 0 M - Bp Ak+ 

0 0 Afc_ -M + Bp 

( 5 ) 

(see Appendix D and E), which assumes the same form 
as the D — 0 BHZ Hamiltonian. Here k± = kx E iky, 
M = Ai 2 — Apoid [see Fig. 1(b)], B = —P/2m*, and A 
characterize the ISOI strength between neighboring mini¬ 
bands with opposite spin. This Hamiltonian obviously 
has a Z 2 topological phase when M < 0, corresponding 
to band inversion. 

Next we employ the eight-band Kane k-p model to cal¬ 
culate the subband structure with SOIs in a 40-nm-thick 


EgS — 



FIG. 2: (color online) (a) Band structure of a 

GaAs/In 2 :Gai_ 2 ;As/GaAs parabolically graded QW from 
eight-band Kane model. The red and blue dashed (solid) 
curves denote the first and second subbands from the 
eight-band Kane model (the reduced four-band model). The 
inset shows the spatial distributions of the first and second 
subbands, (b) and (c): minibands of a parabolically graded 
QW for two different triangular antidot lattices: the lattice 
constants and barrier heights Vo of the antidot lattices are 
a = 17 nm, Vo = 200 meV for (b) and 12.5 nm, 300 meV 
for (c). The minibands are inverted and the minigaps are 
opened. 


GaAs/Ina;Gai_a;As/GaAs PQWs [13, Slj, as plotted in 
Fig. 2(a). The energy difference between the minima of 
the first and second subbands at T point is about 90 meV 
[see Fig. 2(a)]. In order to calculate the miniband struc¬ 
tures caused by an in-plane periodic potential induced by 
the triangular antidot lattice, we reduce the eight-band 
model to an effective four-band k-p Hamiltonian by in¬ 
cluding the lowest 20 electron subbands and 54 highest 
’hole subbands in the QW, to reproduce the energy dis¬ 
persions of the first and second subbands calculated from 
the eight-band Kane model [see Fig. 2(a)]. The param¬ 
eters in the four-band Hamiltonian is given in Appendix 
E. The minibands from the four-band k-p Hamiltonian 
are shown in Figs. 2(b) and 2(c). These minibands orig¬ 
inates from folding the first and second subbands of the 
QW into the first Brillouin zone of the antidot lattice 
[Fig. 1(c)]. By tuning the antidot lattice constant a and 
the potential height Vq, i.e., the etching depth of the an¬ 
tidot lattice, many band inversions appear between these 
minibands, which can be clearly seen in Figs. 2(b) and 
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FIG. 3: (color online) Edge states in nontrivial minigaps of 
a Hall bar structure of the antidot lattice. The Hall bar ori¬ 
entation is along the x axis (a) or y axis (c), as sketched in 
the insets of (b) and (d), respectively. The lower panels (b) 
and (d) amplify the lowest nontrivial minigaps and gapless 
topological edge states. The lattice constant of the antidot 
lattice is a = 17 nm. The blue curves indicate the spatial 
distributions of the edge states. 


2(c). The minigaps between these minibands are opened 
by the ISOI shown in Eq. (g]) [see Figs. 2(b) and 2(c)]. 

To demonstrate that these minigaps are topologically 
nontrivial, we determine the parity of each miniband at 
the four time-reversal invariant momenta [m r, (t = 
0,1, 2, 3) in the first Brillouin zone shown in Fig. 1(c). 
For the lowest N spin-degenerate minibands being oc¬ 
cupied, the Z 2 invariant is given by (—1)*^ = = 

nm=i? 2 m (Tj), where ^ 2 m (Ti) = ±1 is the parity of the 
2TOth occupied miniband at T^. Our calculation gives 
iz = 1 at all the minigaps, which proves the whole system 
is in the quantum spin-Hall phase (see Appendix F). 

Next, we demonstrate the emergence of topological 
edge states upon etching the QW into a Hall bar struc¬ 
ture along two different directions {x axis and y axis). As 
shown in Fig. 3, a pair of topological helical edge states 
appear inside each nontrivial minigap. For example, we 
can see topological helical edge states in the lowest two 
nontrivial minigaps near ^ 186.5 meV and ~ 255 meV, 
respectively. The helical edge state pairs in these mini- 
gaps would lead to higher conductance plateaus as the 
Fermi energy increases by increasing the doping level. 
The helical edge states do not overlap with the bulk 
states, making it possible to be detected experimentally. 

The lowest nontrivial minigaps is quite small (about 
0.5 meV), but the second minigap is larger (about 5 
meV). By tuning the period and potential height of the 
antidot lattice, the nontrivial minigaps can be signifi¬ 
cantly enhanced [see Figs. 4(a) and 4(b)]. For example, 
the lowest minigaps can be enhanced to 5 meV, which is 
already comparable with that in HgTe and InAs/GaSb 


QW systems (~ 10 meV) ii- The second minigaps 
can approach 20 meV, which means the TI phase can 
be realized at liquid nitrigen temperature regime. From 
Figs. 4(a) and 4(b), one can see that the lowest nontriv¬ 
ial minigap is closed as the lattice constant a increases, 
but the second higher nontrivial minigap survives, i.e., 
the TI phase can exist even at large lattice constants, 
e.g., 25 nm. 

Finally, we comment on the experimental detection the 
aforementioned edge states in GaAs/Ina:Gai_a;As/GaAs 
quantum Hall bar (shown schematically in the insets of 
Figs. 3(b) and 3(d)). One way to detect the aforemen¬ 
tioned edge states (shown in Fig. 3) is the standard 
four terminal measurements as demonstrated in previ¬ 
ous works [ 1 , 0 . In contrast to HgTe and InAs/GaSb 
QW systems, there are many pairs of helical edge states 
in our system between these inverted minibands, which 
leads to higher plateaus with increasing the Fermi en¬ 
ergy. Another possible way is microwave impedance mi¬ 
croscopy which makes spatial-resolved nano-scale images 
(< 100 nm) of the conductivity and permittivity of a sam¬ 
ple [ 4 ^. The unoccupied edge states in higher minigaps 
can be detected using the angle-resolved photonemission 
technique Q, which has already been successfully ap¬ 
plied to identify occupied and unoccupied surface states 
in Bi 2 Se 3 and Bi 2 Tea;Se 3 



Lattice Constants(nm) 

FIG. 4: Phase diagrams of the antidot lattice on a 
GaAs/In^^Gai-iAs/GaAs parabolically graded QWs. The 
lowest (a) and second (b) minigap vs. antidot lattice con¬ 
stant a and potential height Vq. Negative (Positive) minigap 
indicates TI phase (normal phase). 
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IV. CONCLUSION 

Our proposal is based on a general analysis about 
the electron orbital motion in TIs. By using the Born- 
Oppenheimer approximation, we find that the fast mo¬ 
tion will induce a spin-dependent gauge field on slow or¬ 
bital motion. Based on this general analysis, we demon¬ 
strate theoretically the TI phase in a conventional 2DEG 
embedded in a symmetric GaAs/Ina;Gai_a;As/GaAs 
PQW, with antidot lattices created by well-developed 
etching technique. The key point is to create a ISOI 
in a symmetric QW, in contrast to conventional SOI in 
asymmetric QWs. This hidden ISOI in symmetric QWs 
induces a spin-dependent effective Lorentz force on the 
electrons, and generates the TI phases in such system. In¬ 
terestingly, such ISOI exists in conventional semiconduc¬ 
tors with a positive bandgap, i.e., normal band structures 
can generate quite large nontrivial gaps approaching 20 
meV. This make it possible to observe the quantum spin 
Hall effect in liquid nitrigen temperature regime. 

So far, all members of TI family are narrow bandgap 
systems containing heavy atoms. Our proposal breaks 
this constraint, and makes it possible to realize TI phase 
in conventional semiconductor 2DEG using the well- 
developed semiconductor fabrication techniques [Tzl - ld^ . 
The presence of the TI phase in PQWs with antidot lat¬ 
tice can largely advance the application of this new quan¬ 
tum state in existing electronics and optoelectronics de¬ 
vices. The general designing principle proposed in this 
work, i.e., the gauge field acting on slow orbital motion 
induced by interband coupling, paves a new way for gen¬ 
erating nontrivial topological phases, such as quantum 
spin Hall phase and even quantum anomalous Hall phases 
by doping magnetic ions, in conventional semiconductor 
2DEGs, and suggests a promising approach to integrate 
it in well developed semiconductor electronic devices. 
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VI. APPENDIX A: EFFECTIVE GAUGE FIELD 
IN SPIN-ORBIT GOUPLED SYSTEMS 


In the spin-orbit coupled system, adopting the Born- 
Oppenheimer approximation, the total Hamiltonian can 
be divided into two parts: 

IL(k, f, a-) = Horb(k, f) -h ffs-o(k, 0-), 


where H^orb(k, f) stands for the intra-band (slow) orbital 
motion part, and i?s-o(k, <t) is the inter-band (fast) spin- 
orbit part. For a given eigenvalue k of the momentum 


operator k, the eigenstates of the spin-orbit part is de¬ 
noted by |n(k)) (n = 1, 2, • • •) and the corresponding 
eigenenergies are e„ (k). We work in the momentum rep¬ 
resentation of the orbital part and expand the eigenstate 
of H(k,r,&) in this representation, |'I'(k)) = (k|'I'), as 

l^(k)) = ^(A„(k)|x„(k)). 

n,k 


In the momentum representation, we have f = iVk and 
k = k. Substituting into iL(k, f, 6-)|4>(k)) = i?|4'(k)), 
we have 

^iJm,«(k)(/)„(k) = 


where 

-^m,n(k) = ^ 

m,n^m (k) -H (x™(k)|iJorb(k,f)|x„(k)) 

— ^m,nCm(k) 4- .Horb(k, f A^ yj(k)) 

contains a pure gauge A,.,^_„(k) = -i(x,.,^(k)|Vk|x„(k)). 
By now the above equation is still exact. Now we make 
the Born-Oppenheimer approximation and consider adi¬ 
abatic transport, i.e., neglect the off-diagonal coupling 
between different spin-orbit energy bands, to arrive at 
the single-band description 

^«(k,f'„)<()„(k) = E(l)n{]<.), 

where the effective single-band Hamiltonian on the or¬ 
bital motion 

HniKK) = en(k) -4iLorb(k,f'„), 

where = f„ — A„(k) contains an effective gauge field 
for the slow orbital motion: 


An(k) = -i (Xn(k)l Vk Ix«(k)) . 


Specifically, for the BHZ model of 2D TIs in the pres¬ 
ence of an in-plane uniform electric field, the slow orbital 
part is 


iLorb(k, r) = C — Dk'^ — eE • r, 
and the fast spin-orbital part iLs-o(k, <t) is 


f M - Bk^ Ak+ 0 0 

Ak- -M -4 0 0 

0 0 M - Bk"^ Ak_ 

Vo 0 Ak+ -M + Bk"^ 


and 


A - s -A. 
^x,n — *’^ 2^2 


1 + In 


Bk^ - M 


-h (M- Hfc2)2]l/2_ 



ZK 


Bk^ -M 


[A2/c2 -h (M- Hfc2)2]l/2 
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where s„ = ±1 for spin-up/down block, and = ±1 for 
the electron in conduction/valence band (n = 1 , 2 ,3,4). 
The effective vector potential leads to the non-trivial ef¬ 
fective gauge field with the strength 


-^a:y,n(k) — i\x ,y ] — -^y\ 


= (Vk X A)^ = A„ 


A^{M + Bk‘^) 


2[A2A:2 + (M-Bfc2)2]3/2^ 


where A„ = x = ±1 (n = 1,2,3,4). Within the 
Born-Oppenheimer approximation, the equation of mo¬ 
tion for the n-th band can be written as 

— TTT; T Fxy,n{k^ky^ 


Vn = 


hdkx 

dH„ 


hdk 

ki = eEi/h 


Exy,n{k^kx ^ 


we can see that the gauge field strength F^y^n = (Vk x 
A) 2 acts as a Lorentz force in the /c-space, acting on spin- 
up and spin-down electrons in opposite directions, which 
is perpendicular to electron momentum. 


VII. APPENDIX B: EFFECTIVE SPIN-ORBIT 
COUPLING IN A QUANTUM WELL 


For a symmetric QW grown along (001) direction (the 
z axis), effective spin-orbit coupling exists between sub¬ 
bands with opposite parities. This effective spin-orbit 
coupling comes from interband coupling and can be un¬ 
derstand by reducing the 8 x 8 Kane Hamiltonian to a 
2x2 effective Hamiltonian [50j |. 

To the first order of fc, the 8 x 8 Kane Hamilto¬ 
nian in the basis (i 15,t), * l-S”,4-), |3/2,1/2), |3/2,—1/2), 
13/2,3/2), |3/2,-3/2), |l/2,l/2), |l/2,-l/2)) around 
the r point is 


-ffsxS — 


He Hex\ 

Hi Hx ) 


where He = ecl 2 x 2 and Hy = Cyhxi ® £ 5 ^ 2 x 2 are 2 x 2 
and 6 x 6 diagonal part for conduction and valence bands, 
and the 2 x 6 matrix 


y/2Pk^. 

-Pk^ 


n 

Pkr. 

-Pk^ 

Ve 

Ve 

V2 


Ve 

Ve 

Pk+ 

V2Pkx 

n 

Pk- 


Pk^ 

Ve 

Ve 


V2 

vA 

Ve 



represents the interband coupling. Specifically, Ci = 
h^k^ / {2m) + Vi is the kinetic energy plus the total po¬ 
tential for the conduction/valence/spin-split (z = c/v/s) 
bands, with Ve — Vy = Eg the band gap and Vy — 
Vs = Aq the band off set. k± = kx i iky and P = 
—i{h/mQ) {Se\px\Xv) parameterize the interband cou¬ 
pling. 

The eigenvalue problem can be expressed as 


( He He, 

\Hl Hy 




where pe is a two-component spinor for conduction 
bands and is a six-component spinor for valence 
bands. Since we focus on the conduction bands, 
^Py = {e - Hy)~^ Hl'tj)e can be eliminated and gives 
the effective Schrodinger-type equation Elgsipe = Sipe, 
with Hes = He + Hey {e — Hy)~^ Hi for conduction 
bands. Without loss of generality, we assume the 
QW is non-uniform only along the 2 ; direction, e.g., 
a PQW. By straightforward algebra, we have H^g = 

= Ei=y,si^-mPl)[ie-e,iz))-^P{z),kl, 

where i/Jg and H^g represent the effective spin-orbit cou¬ 
pling between the spin up and down electron. 

Since we focus on the lowest conduction subbands, 
we have e — ey{z) ~ Eg{z) and e — es{z) ~ Eg{z) -I- 
Ao(z) = Egi{z). Because Eg{z) and Egi{z) are much 
larger than the subband energies in the wide QWs under 
consideration, we keep the zero-th order terms E~^{z) 
and Egi{z)~^ in the expansion, and project the spin- 
orbit coupling operator into the two lowest spin- 
degenerate subbands (|xi(-z))j \X 2 {z))) to obtain the ISOI 
V'''x{kxCk^ + kya^), where rj is given in Eq. ([4]), cr* denotes 
the real electron spin, and refers to the Pauli matrix 
describing the the subband index. 


VIII. APPENDIX C: BAND EDGE WAVE 
FUNCTIONS IN FOLDED BRILLOUIN ZONE 


We consider a PQW in the presence of an antidot lat¬ 
tice, which can be generally described by a potential 
V{r) = X(q V(q)e*'i''’ with the lattice periodicity. 

For a triangular antidot lattice, the reciprocal lat¬ 
tice vectors in the hexagonal Brillouin zone are Qi = 
(2V3a)(3,0), Q 2 = (27r/3a) (-3/2,-3^3/2), Q 3 = 
— (Qi + Q2)- The envelope functions of the low¬ 
est miniband at the band edge (fc = 0, F point) is 

Mr,i(i’) = 1 - SLi ‘^'rnV{Qi)/ffQ'^ cos(Q/ • r). For 


higher minibands, their envelope functions ur,rt(r) {n = 
2,3,4, 5,6, 7) at the band edge are linear combinations 
of the six wave vector components 

e.g., ■ur. 2 (r) oc sin(Qi • r) and ■ur. 7 (r) oc 
cos(Qi • r). The most important minibands are 
Uk,i{r), Uk, 2 {r) and Uk, 4 {r): the lowest nontrivial mini¬ 
gap occurs between Uk,i{r) and Uk,2{f'), and the second 
nontrivial minigap occurs between UkpD and Uk, 4 {r). 


IX. APPENDIX D: EFFECTIVE BHZ 
HAMILTONIAN NEAR THE V POINT 

The lowest two subbands XiD and X 2 {z) in a PQW 
have even and odd opposite parities, an effective spin- 
orbit interaction {xi{z)\ri{z){kxcr^ + kya^)\x 2 {z)) ap¬ 
pears. When the Brillouin zone is folded by the tri¬ 
angular anti-dot lattice, the lowest nontrivial minigap 
appears between the miniband pair \xi{z),Uk, 2 i'r)) and 
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\X 2 {z),Uk,i(r)), i..e., the second miniband of the first sub¬ 
band and the first miniband of the second subband. The 
second nontrivial minigap appears between the miniband 
pair 1 x 2 ( 2 ),■Ufc, 2 (?')) and \xi{z),Uk,i{r)), i.e., the second 
miniband of the second subband and the fourth mini¬ 
band of the first subband. To obtain an effective Hamil¬ 
tonian near each minigap, we project the Hamiltonian 
H = /2m-\- IS.i 2 Tz + r]Tx{kx(J^ + kyCj'^) -\- V(x, y) onto 

the corresponding miniband pair and obtain an effective 
BHZ model Eq. © in the basis |^+,t), |'0-,4.), |V^+,4.), 
where \4’+) is the miniband above |'0_) by 2M 
at the r point, B = —characters the band dis¬ 
persions with the effective mass m* near the band edge, 
and A characterize the intersubband spin-orbit coupling. 
At the r point 

(V'-hI v{z){k^a-‘ + kya^) \ip-) 

= (X1(2)(2)| |X2(1)(2)) 

■ {uk,2{2){r) \ - idyUy) |wfe,i(4)(?')) 

= A {kxO-^ + kya^). 

The accurate coupling strength can be estimated by nu¬ 
merical calculating based on the eight-band Kane model. 

For BHZ model, a Z 2 topological transition from the 
normal phase to the topological insulator phase would 
occur when M = A 12 — Apoid [see Fig. 1(b)] changes 
sign from positive to negative, which can be controlled 
by adjusting the lattice constants and etching depths of 
antidots. 


X. APPENDIX E: THE EFFECTIVE 
HAMILTONIAN REDUCED NUMERICALLY 
FROM THE EIGHT-BAND KANE MODEL 


At the r point, the wave functions in the eight-band 
k • p Hamiltonian are 


/ eI"^^ 

(^)\ 


(^) 

p(m) 

iz) 


(z) 

V Fi’”) 

iz)/ 


which can be obtained by solving the secular equation 

Considering the two lowest electron subbands, we ob¬ 
tain the effective two-dimensional Hamiltonian by aver¬ 
aging the 2 : component in the Hamiltonian 

i7eff(k||) = (vI/(z)|7J|d/(z)), 


where the matrix element of the Hamiltonian is 

= ( 2 )). 

i,I=l 

The Hamiltonian can be divided into 
H = + H' 


H 

Then we have 


' (q, fcz) = a (q) + ^ (q) ^2+7 (q) kl. 




E 

*j'=i 




+ iz)\[Piq)],^k\F^^^ ( 2 )) 

+ {Ft^ {z)\bi<D],,kl\Fj-^ (z)). 


The contribution of the subbands other than the two 
lowest electron subbands should also be considered in the 
reducing process, which can be done by using Lowdin 
perturbation theory. We include the lowest 20 electron 
subbands and 54 highest hole subbands in the QW re¬ 
spectively and divide them into the weakly coupled sub¬ 
sets and 5*2. The set S*! includes the two lowest elec¬ 
tron subbands |xi) and |x 2 ), the other subbands are in¬ 
cluded in the set S 2 - The Hamiltonian is reduced into 
set Si using the Lowdin perturbation method, 


^(2) _ \ ' jzf' rr' 


+ 


Em — El Em' — El 


where the indices m correspond to states in the set A, 
the indices I correspond to states in the set B, and 

Hmi = (V’mj-ff' \ipi) ■ 


Finally we obtain the effective two-dimensional Hamil¬ 
tonian in the basis |xi,t), 1X2,4-), Ixi,)-), 1X2, t): 


/ El + Bi ■ k^^ Afc+ 0 

Ak_ E2 + B2-kl 0 
0 0 Ei+Bi-k^ 

Vo 0 Ak+ 


0 

0 

Ak- 

E 2 -l- B 2 • fcjj 


where 


El = 0.05351eV, 

Bi = 0.86571eV • nmV 
E 2 = 0.14396eV, 
i ?2 = 0.67969eV • nm^, 
A = 0.01041eV • nm. 













In order to examine the validity of the 4-band 
Hamiltonian, we plot the band structure of the 
GaAs/Ina;Gai_a,As/GaAs PQW calculated by the 4- 
band model and compare it with the eight-band k • p 
model (see Fig. 2a in the manuscript). One can see 
clearly that the band structure obtained from the 4- 
band model [the solid lines in Fig. 2a] is in good agree¬ 
ment with that obtained from the eight-band model [the 
dashed lines in Fig. 2a]. 

XI. APPENDIX F: VERIFICATION OF 
NON-TRIVIAL Za TOPOLOGICAL INVARIANT 

Topological insulators with dissipationless edge states 
and ordinary insulators are distinguished by different Z 2 
invariants. For 2D systems, Fu and Kane [ll| have shown 
that the Z 2 invariant can be determined from the parity 
of the occupied band at the four time-reversal invariant 
momenta in the Brillouin zone. The Z 2 invariant ^ = 0,1, 
which distinguishes the quantum spin-Hall phase in two 
dimensions, is given by 

(-ir = n.<5., 

~ nm=1^2m (Fj) , 


where ^ 2 m (Fi) = ±1 is the parity eigenvalue of the 2mth 
occupied energy band at the time-reversal invariant point 
Fi, which shares the same eigenvalue ^ 2 m = ^ 2 m-i with 
its Kramer degenerate partner. The four time-reversal 
invariant points = {n-iQi + n 2 Q 2 )/2, where 

'ni,n 2 = 0,1. The calculated parity eigenvalue of the 
2mth (to = 1, 2, 3,4, 5,6 ) occupied energy band at F^ are 
listed: 


ra 

^2m (hoo) ^2rr 

1 (roi) ^2rr 

(^ 10 ) ^2rr 

(rii) 

1 

+ 

- 

- 

- 

2 

0 

+ 

-h 

+ 

-h 

4 

+ 

+ 

+ 

-h 

5 

+ 

-h 

-h 

-h 

6 

- 

- 

- 

- 


From the above calculation, we can confirm that the 
Z 2 invariant v = 1 a.t the 2TOth (to = 2, 5, 6 ) occupied 
band where the minigaps open, and the system enters 
the TI phase and the dissipationless edge states appear. 
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